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Metastable  Molecules  and  Other  Energehc  Structures: 
T HE  Contribution  of  Theory 


Rodney  J.  Bartlett 
Principal  Investigator 
Quantum  Theory  Project 
University  of  Florida 
Gainesville  FL  32611 


I.  Introduction 

Our  prior  work  for  AFOSR  has  focused  on  developing  electron-correlated 
ab  initio  \heor^,  implementing  new  computational  methods,  and  predicting  new, 
high-energy  molecules.  Our  recent  (last  three  years)  achievements  in  this  field 
are  summarized  in  the  published  papers  and  conference  presentations  listed  in 
Sections  IV  and  V.  We  have  emphasized  a  rational  theoretical  design  of  new 
energetic  metastable  molecules  that  is  based  on  highly  accurate  electron- 
correlated  coupled-cluster  theory  and  directed  toward  a  search  of  "synthesizable" 
species  that  have  a  high  heat  of  formation  (large  specific  impulse)  and  can  be 
reasonably  stable. 

We  have  Investigated  the  molecular  mechanisms  of  the  initial  steps  in 
exothermic  decomposition  (initial  steps  of  ignition  and  explosion),  as  well  as 
possible  synthetic  roots  in  the  synthesis  of  the  desired  metastable  isomers.  The 
barriers  and  the  transition  structures  of  these  reactions  provide  clear  clues  to 
possible  synthetic  approaches  and  to  the  stability  of  all  final  and  intermediate 
species,  as  well  as  provide  the  spectroscopic  "fingerprints"  in  assessing  the 
pathways  the  actual  reactions  will  follow. 

Based  on  the  goal  of  developing  theoretical  methods  that  can  provide 
sufficiently  high  accuracy  to  be  competitive  with  what  experimental  methods  can 
achieve  for  the  characterization  of  molecular  structure  and  properties,  in  this 
work  we  have  formulated  and  implemented  in  the  ACES  II  program  the  following 
new  "tools": 

♦  ROHF-CC  and  ROHF-MBPT  methods  for  the  treatment  of  open-shell 
systems  that  might  be  subject  to  spin  contamination  in  UHF 

♦  Analytical  gradients  for  ROHF-CC  and  ROHF-MBPT 

♦  The  equation-of-motion  coupled-cluster  (EOM-CC)  approach  for 
excited,  electron-attached,  and  ionized  states 

♦  Triple  excitations  for  excited-state  energies 
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♦  Simplified  (partitioned)  EOM-CC  methods  for  application  to  large 
molecules 

♦  EOM-CC  method  for  computations  of  properties,  including  NMR 
coupling  constants 

The  synthesis  of  new  high-energy  density  materials,  rocket  propellants, 
and  explosives  has  reached  a  limit  where  the  design  of  a  new  HEDM  will  require 
a  detailed  knowledge  of  molecular  mechanisms  and  eventually  of  rate  constants 
of  reactions  that  lead  to  the  formation  or  to  the  decomposition  of  metastable 
molecules.  In  general,  the  increase  in  energy  content  (heat  of  formation)  would 
rather  decrease  the  stability  of  the  substances,  if  there  are  no  particular 
electronic  mechanisms  (e.g.,  aromatic  stabilization  or  other  delocalization  effects 
arising  from  symmetry)  or  geometrical  factors  (e.g.,  separation  of  reactive 
groups  by  bulky  substituents  or  by  solvent  molecules),  which  stabilize  the 
energetic  molecules.  Nevertheless,  one  should  remember  that  approaching  the 
limit  of  chemical  energy  storage  will  also  decrease  the  effect  of  electronic 
stabilization.  This  pertains  not  only  to  the  final  energetic  molecules,  but  also  to 
the  reaction  intermediates,  which  may  cause  the  decomposition  of  the  final 
product  at  the  moment  of  its  formation.  Thus  only  a  detailed  knowledge  of  the 
reaction  mechanisms  can  help  to  synthesize  and  to  find  the  requisite  conditions 
to  prevent  immediate  decomposition  of  very  energetic  species. 

II.  Summary  OF  Principal  Results 

Ab  Initio  computations  of  the  heats  of  formations  and  structures  of  known 
or  easily  predicted  unknown  molecules  have  become  routine  using  programs  like 
ACES  II,  Gaussian  or  GAMESS.  ACES  is  a  user-friendly  program  which  allows  it  to 
be  employed  as  a  tool  by  an  experimental  chemist  without  requiring  a  special 
knowledge  of  theory,  yet  it  incorporates  the  newest  developments  we  are 
making  in  CC/MBPT  and  their  excited-state  extensions.  This  has  enabled  us  to 
make  comprehensive  studies  of  new  energetic  molecules  including  possible 
synthetic  routes  and  decomposition  reactions.  Besides  incorporating  our  new 
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methods,  we  are  also  trying  to  improve  the  performance  of  the  ACES  II  program 
with  the  intention  of  providing  it  as  an  easily  used  universal  tool  for  structural, 
energetic,  and  spectroscopic  studies  of  metastable  molecules. 

A.  New,  Potential  Monopropellants  and  Oxidizers 

One  of  our  objectives  has  been  to  pursue  systematic  studies  of  small 
energetic  molecules  that  exhaust  the  isomeric  possibilities.  In  this  way,  we  have 
explored  the  CN2O2  [25]  and  CN2O3  [31]  systems,  while  also  considering 
selective  dimers  and  trimers  [26]. 

1.  Nitryi  cyanide,  NC-NO2  {1}  [25]  is  composed  of  two  energetic  groups, 
CN  and  NO2.  The  estimated  large  values  of  the  heat  of  formation 
(AH°3oo  =  60  kcal  /  mol)  and  of  the  decomposition  energy  of  1 
(1  Nj  +  CO2;  AE  =  150  kcal  /  mol)  make -this  species  potentially  interesting  as 
a  high-energy  molecule.  Nittyl  cyanide  can  be  a  reasonably  stable  molecule  as  its 
C-N  bond  dissociation  energy  (59  kcal/mol)  and  the  barrier  to  decomposition  into 
N2  and  CO2  (54  kcal/mol)  are  rather  high,  being  comparable  to  nitromethane. 
Nitrosoisocyanate,  ONNCO  {2},  earlier  studied  as  a  collision  complex  in  the 
reaction  of  NCO  and  NO  (M.C.  Lin  and  C.  F.  Melius,  J.  Phys.  Chem.  1993,  97, 
9124)  is  the  most  energetically  favorable  CN2O2  isomer,  but  its  18  kcal/mol 
unimolecular  dissociation  barrier  is  very  low.  Thus  2  can  only  be  observed  as  a 
short-lived  intermediate.  Nitrosofulminate,  ONCNO  {3},  and  its  bicyclic  isomer 
■  {4}  are  even  more  energetic  than  1,  but  they  are  apparently  less  stable. 
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2.  2,4,6-trinitro-l,3,5-triazine  (TNTA)  {5}  [26]  is  the  cyclic  trimer  of  1. 
The  formation  of  molecular  complexes  (e.g.,  di-,  tri-,  oligo-,  and  polymers)  is  a 
conceivable  way  to  increase  density  and  stability  and  to  improve  the  material 
properties  of  propellants  and  explosives.  Trinitrotriazine  has  an  intriguing  history. 
As  suggested  in  1907  (Finger,  H.,  J.  Prakt.  Chem.  75,  103  [1907]),  one  might 
expect  the  synthesis  of  TNTA  from  industrially  available  cyanuric  chloride  (2,4,6- 
trichlorotriazine)  and  AgN02,  however,  this  does  not  occur.  All  previous  efforts  to 
synthesize  1  failed,  and  no  further  information  about  preparation  or 
characterization  of  TNTA  is  available  in  the  literature.  Our  computational  results 
show  that  TNTA  can  be  reasonably  stable,  if  the  proper  synthetic  route  is  found. 


The  minimum  energy  C3  structure  of  TNTA  computed  at  the  MBPT(2)/6- 
31G*  level  has  a  planar  aromatic  six-membered  ring  with  the  NO2  group  rotated 
out  of  the  ring  plane  by  40°.  The  out-of-plane  rotation  of  NO2  groups 
compromises  the  conjugation  of  the  nitro  groups  and  the  aromatic  six-electron 
ring  system  (maximum  for  the  planar  form)  and  repulsion  of  negatively  charged 
oxygens  and  ring  nitrogens  (minimum  for  the  90°  rotation).  The  46  kcal/mol 
heat  of  formation  of  TNTA  can  be  obtained  from  an  isodesmic  reaction,  and  the 
estimated  heat  of  formation  of  TNTA  in  the  solid  state  is  24  kcal/mol  lower  than 
its  gas  phase  value.  Once  we  have  an  estimate  of  AH°f  for  TNTA,  its  304 
kcal/mol  energy  release  in  decomposition  into  N2  and  CO2  (AEr),  specific 
decomposition  energy  (1410  cal/g;  AEr  divided  by  the  molecular  mass)  and 
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specific  impulse  (269  s  in  the  atmosphere  and  294  s  in  the  vacuum)  can  be  easily 
derived.  Our  estimates  of  AEr  and  Isp  for  TNTA  demonstrate  its  performance  to 
be  similar  to  that  of  RDX  (hexahydro-l,3,5-trinitro-l,3,5-triazine),  but  the 
estimated  density  of  TNTA  (2.1  g/cm^)  is  higher  than  in  RDX  (1.82  g/cm^). 

3.  Nitrofulminate,  0NC-N02  {6}  [31],  N-oxide  of  nitryl  cyanide,  promises 
to  be  a  very  powerful  oxidizer,  if  a  method  for  its  synthesis  is  found.  Regarding 
its  high  exothermicity  in  decomposition,  nitrofulminate  is  suggested  as  a 
potential  energetic  oxidizer.  The  estimated  gaseous  heat  of  formation  of 
nitrofulminate  is  71  kcal/mol,  and  it  has  very  high  exothermicity  in 
decomposition: 

O2N  -  CNO  ^  CO2  +  N2  + 1  /  2O2;  AE  =  -165  kcal  /  mol. 

Another  CN2O3  isomer,  nitroisocyanate,  OCN-NO2  {7}  is  a  collision 
complex  in  the  gas  phase  reaction  of  NO2  and  NCO  radical.  The  reaction  of  NO 
and  NCO  has  a  similar  mechanism  that  involves  the  formation  of  the  unstable 
nitrosoisocyanate,  OCN-NO.  Among  other  reasons,  an  interest  in  the 
experimental  and  theoretical  studies  of  reactions  of  NCO  with  NOx  occurs 
because  of  their  importance  for  the  process  of  effective  removal  of  NOx  from 
exhaust  gases— the  rapid  reduction  of  nitrogen  oxide  (RAPRENOx)  (Perry,  R.  A., 
Siebers,  D.  L,  Nature  324,  657  [1986]). 


4.  Trinitrotriazine  trioxide  (TNTATO)  {8}  [31]  is  an  analog  of  triazine,  but 
the  increased  number  of  oxygen  atoms  offers  7  as  a  potential  oxidizer  with  a 
high  heat  of  formation  (AH°f(g)  =  108  kcal/mol).  As  estimated  in  composite 
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liquid  and  solid  propellants,  the  performance  of  TNTATO  is  comparable  to  that  of 
the  most  energetic,  experimentally  known  oxidizers.  The  formation  of  two  other 
cyclic  triazine-type,  (CN203)3  trimers  {9, 10}  which  have  lower  energy,  may 
complicate  the  synthesis  of  8  or  decrease  its  stability. 


5.  The  20,  eight-electron-four-atom  metastable  molecules  represent 
another  type  of  grouping  that  might  be  explored  to  assess  some  characteristics 
of  energetic  species.  This  was  suggested  as  an  extension  of  our  previous 
research  on  the  nitrogen  dimer  [21].  We  have  found  that  the  alternative 
decomposition  pathway  of  tetraazatetrahedrane  {11}  via  the  planar  four- 
membered  anti-aromatic  tetrazacyclobutadiene  {12}  has  a  higher  barrier  than 
the  direct  decomposition  into  two  N2  molecules. 

A  few  other  molecules  of  the  series  have  interesting  structural,  energetic, 
and  spectroscopic  features  that  can  be  of  potential  interest  for  further 
experimental  study  and  possible  applications.  Thus  two  CN2O  isomers, 
diazirinone  {13}  and  nitrosyl  cyanide  {14},  are  high-energy  molecules.  Although 
14  has  a  somewhat  higher  energy,  it  is  known  experimentally;  13  is  unknown. 
Other  unknown  molecules  with  intriguing  properties  are  fluorodiazaboririne  {15}, 
ethylenedione  {16},  difluorodiborene  {17},  and  nitrosoboroxide  {18}. 
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B.  Methodological  Advances 

We  have  made  substantial  progress  in  the  following  areas: 

♦  Investigation  of  the  equation-of-motion  coupled-cluster  approach  for 
excited,  electron-attached,  and  ionized  states; 

♦  Introduction  and  assessment  of  the  effects  of  tripie  excitations  on 
excitation  energies; 

♦  Study  of  simplified  (partitioned)  EOM-CC  methods  for  appiication  to  large 
molecules; 

♦  Consideration  of  the  application  of  EOM-CC  to  properties,  including  NMR 
coupling  constants. 

1.  The  electronic  excitation  (EE)  EOM-CC  limited  to  single  and  double 
excitations  (i.e.,  EOM-CCSD)  has  now  been  applied  to  a  range  of  difficult 
systems,  including  the  molecules  CI-O-CI  [14],  formaldehyde,  acetylaldehyde, 
and  acetone  [12],  the  polyenes  (ethylene,  butadiene,  clyclopentadiene)  [24],  and 
benzene  and  the  azabenzenes  [30].  It  is  clear  from  these  studies  that  with  well 
polarized  and  diffuse  basis  sets,  EE-EOM-CCSD  is  accurate  to  within  about  0.3 
eV,  except  in  cases  where  the  excited  state  is  dominated  by  a  double  excitation. 
Whereas  the  errors  are  small  for  Rydberg  states,  they  tend  to  also  be  large  for 
some  valence  excited  states.  The  above  points  to  a  need  for  triple  excitations  to 
rectify  the  problem  with  doubly  excited  states  and  with  certain  valence  states. 

2.  A  hierarchy  of  triple-excitation-corrected  EOM-CC  methods  has  been 
introduced  in  the  last  two  years  of  this  project  [9,  23].  There  are  iterative 
approximations  based  on  our  CCSDT-1  [Y.  S.  Lee  and  R.  J.  Bartlett,  J.  Chem. 
Phys.  134,  126  (1984)]  and  CCSDT-3  [J.  Noga,  R.  J.  Bartlett,  and  M.  Urban, 
Chem.  Phys.  Lett.  134,  126  (1987)]  methods,  as  well  as  their  corresponding  non¬ 
iterative  approximations,  EOM-CCSD(T)  [9]  and  EOM-CCSD(T)  [23].  The  non¬ 
iterative  methods  depend  on  two  and  three  n^N"^  steps,  respectively,  rather  than 
n^N\Nit  as  In  the  interative  methods,  so  they  can  be  applied  to  larger  molecules. 
The  EOM-CCSDT-1  and  (T)  methods  have  been  shown  to  offer  much  improved 
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results  for  double  excited  states,  while  they  do  little  for  the  valence  excited 
states.  However,  the  EOM-CCSDT-3  offers  substantial  improvement  for  both, 
showing  exceptional  agreement  in  comparisons  to  the  full  Cl  limit  for  several 
small  molecules  [23].  Most  of  this  accuracy  is  retained  by  the  much  less 
expensive  EOM-CCSD(T).  For  the  polyenes  [24]  and  azabenzenes  [30],  the 
average  error  for  the  latter,  compared  to  experimentally  established,  vertical 
excitation  energies,  is  about  1.0  eV.  It  does  overestimate  the  triple  correction 
compared  to  EOM-CCSDT-3  for  some  difficult  cases,  however. 

3.  For  large-scale  applications  even  EOM-CCSD  is  highly  demanding,  as 
nonsymmetric  matrices  of  rank  exceeding  10^  must  be  partially  diagonalized 
iteratively,  to  extract  the  eigenvectors  and  eigenvalues.  Consequently,  we  have 
investigated  approximations  based  upon  partitioning  such  large  matrices  to 
greatly  reduce  the  computational  complexity  while  retaining  most  of  the  accuracy 
of  the  full  EOM-CCSD.  In  particular,  we  introduced  the  p-EOM-CCSD  and  its 
second-order  perturbation  theory  analogue,  p-EOM-MBPT(2)  [17].  The  latter 
requires  only  an  n^N^  algorithm  compared  to  n^N"*  for  the  EOM-CCSD,  and  serves 
as  the  analogue  of  the  highly  successful  MBPT(2)  treatment  for  ground  states.  In 
applications  to  ethylene,  butadiene,  formaldehyde,  and  acetylaldehyde,  the  p- 
EOM-MBPT(2)  method  shows  an  average  error  of  0.12eV  compared  to  the  full 
EOM-CCSD  results.  However,  for  significantly  large  systems,  such  a  calculation  is 
one  to  two  orders  of  magnitude  faster  while  also  being  far  less  demanding  of 
computer  resources. 

4.  The  electron-attached  (EA)  EOM-CC  was  formulated  and  implemented 
recently  [10].  Besides  the  obvious  application  to  electron  affinities,  it  is  ideally 
suited  to  describing  systems  that  differ  by  one  electron  from  a  closed  shell.  An 
example  is  offered  by  our  work  on  the  Li  embedded  cryogenic  cluster  problem. 
For  such  a  system  (consisting  of  an  alkali  metal  atom,  M,  plus  a  closed  shell 
species— Hz,  Nz,  He,  Ne,  Ar,  etc.),  we  can  study  the  excited  states  by  first 
performing  CC  calculations  on  the  M"^  -i-  Hz  system  and  then  using  EA-EOM-CC  to 
systematically  add  an  electron  into  each  of  the  unoccupied  orbitals  in  the  cation. 
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This  conceptually  appealing  route  also  has  numerical  advantages,  as  there  are 
fewer  coefficients  to  determine  than  in  EE-EOM-CC,  yet  it  provides  a  direct,  fully 
correlated  calculation  of  the  comparatively  small  energy  difference.  Such  an 
approach  benefits  from  the  large  similarity  between  the  correlation  problem  for 
the  cation  and  the  neutral  molecule.  For  the  Li  +  H2  system,  our  answers  are 
exceptional,  although  the  simultaneous  description  of  the  weak  ground  state 
Van-der-Waals  interaction,  in  addition  to  that  for  the  excited  states,  makes  the 
problem  exceedingly  difficult.  The  EA-EOM  problem  is  also  pertinent  to  our  newly 
proposed  STEOM-CC. 

5.  A  very  convenient  treatment  of  second-order  properties  is  a  byproduct 
of  the  EOM-CC  approach.  Those  of  primary  interest  here  are  the  spin-spin 
coupling  constants  since,  unlike  NMR  chemical  shifts,  they  have  not  been 
adequately  described  in  prior  theory.  Furthermore,  they  offer  characteristic 
values  that  can  identify  various  bonding  situations  likely  to  be  encountered  in 
atom-doped  cryogenic  solids.  To  enable  theory  to  handle  this  problem,  we 
introduced  new  EOM-CC  property  methods  [3,  4]  in  1994.  This  work 
demonstrates  that  it  is  necessary  to  include  all  four  of  the  contributions  to  the 
coupling  constant  (i.e.,  the  Fermi  contact),  spin-dipolar,  paramagnetic  spin  orbit, 
and  the  diamagnetic  terms,  if  predictive  accuracy  is  to  be  attained.  Good  basis 
sets  and  high  (CC)  levels  of  electron  correlation  are  required.  After  calibration 
[19],  where  we  determined  average  errors  of  4  Hz  (or  less  than  about  10%  for 
the  coupling  constants),  we  had  a  tool  to  make  "predictive"  applications.  We  did 
so  for  the  fascinating  carbocations,  which  can  assume  non-classical  strucutres. 
First  we  studied  the  ethyl  carbenium  ion  [13],  This  prototypical  primary  alkyl 
cation  has  never  been  observed  experimentally,  but  our  calculations  provide  a 
coupling  constant  of  13  ±  7  Hz  for  the  bridged  H-C,  which  would  serve  as  a 
fingerprint  for  the  non-classical  H,  as  the  usual  H-C  constant  is  about  120  Hz. 
Similarly,  the  bridged  C  to  a  normal  C  coupling  constants  in  the  2-norbornyl 
cation  (the  system  studied  extensively  by  G.  Olah)  have  never  been  determined 
experimentally.  Our  calculation  [22],  which  provides  excellent  agreement  with 
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the  constants  that  were  determined,  predicts  -5  Hz  compared  to  the  normal  C-C 

constants,  which  are  about  +33  Hz. 
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Spectroscopy:  Vibrational,  Electronic,  and  NMR,"  Korean  Chemical  Society, 
Pusan,  Korea. 

(16)  September  1995— "Equation-of-Motion  Coupled-Cluster  Method  for  Excited, 
Ionized,  and  Electron-Attached  States,"  Molecular  Quantum  Mechanics: 
Methods  and  Applications  Conference,  University  of  Cambridge,  Cambridge, 
England. 

(17)  August  1995— "Excited  States  Using  the  Equation-of-Motion  Coupled-Cluster 
Method,"  12‘^  Canadian  Symposium  on  Theoretical  Chemistry,  University  of 
New  Brunswick,  Fredericton,  Canada. 

(18)  June  1995— "Analogues  of  Tetrahedral  N4  as  Potential  High-Energy 
Molecules,"  Air  Force  High  Density  Materials  Contractors'  Conference, 
Woods  Hole,  MA. 
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(19)  May  1995— "Equation-of-Motion  Coupled-Cluster  Methods  for  Excited, 
Ionized,  and  Electron-Attached  States,"  Pople  Symposium,  Northwestern 
University,  Evanston,  IL. 

(20)  May  1995— "Equation-of-Motion  Coupled-Cluster  Methods  with  Application 
to  BCI3,"  Air  Force  Office  of  Scientific  Research  Contractors'  Meeting, 
Wright-Patterson  Air  Force  Base,  CA. 

(21)  January  1995— Air  Force  Office  of  Scientific  Research  HEDM  Theory- 
Synthesis  Workshop,  Palmdale,  CA. 

University  Colioquia  Presented 

February  1997— Tulane  University,  New  Orleans,  Louisana 
December  1996— University  of  Torino,  Torino,  Italy 
April  1996— University  of  North  Texas,  Denton,  TX 
April  1996— Southern  Methodist  University,  Dallas,  TX 
April  1996— University  of  Texas  at  Arlington,  Arlington,  TX 
October  1995— Lehigh  University,  Bethlehem,  PA 
March  1995— Iowa  State  University,  Ames,  lA 

V.  Honors,  Awards 

David  L.  Williams  Term  Professorship 

Dr.  Rodney  J.  Bartlett  has  been  designated  the  1997-1998  David  L. 
Williams  Term  Professor  in  the  College  of  Liberal  Arts  and  Sciences  at  the 
University  of  Florida.  The  award  is  given  in  recognition  of  outstanding  scholarship 
and  teaching  and  carries  with  it  a  one-year  supplement  of  $5,000  in  salary  and 
$1,000  in  research  support. 

University  of  Florida  Research  Foundation  Professorship  Award 

Dr.  Rodney  J.  Bartlett  has  been  awarded  a  1998  University  of  Florida 
Research  Foundation  Professorship.  These  are  three-year  professorships  to 
recognize  faculty  who  have  distinguished  records  of  research  and  scholarship 
that  clearly  will  be  sustained.  The  award  includes  a  salary  bonus  and  research 
support  money. 
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